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THE RESPONSE OF THE THERMOSPHERIC DENSITY 
TO AURORAL HEATING DURING GEOMAGNETIC 

DISTURBANCE 

H. Volland and H. G, Mayr 
ABSTRACT 

Assuming an impulse type heat input into a small band of latitude within the au-> 
rora ovals during local night; which shall simulate a heat input during geomagnetic 
disturbances, the corresponding response of the thermospheric density has been 
calculated. The result in terms of a series of spherical harmonics shows that the 
components with large wave domain numbers (n, m) decay rapidly within the first 
hour after the onset of the geomagnetic storm while the two zonal components 
(0, 0) and (2, 0) and the two associated components (1, 1) and (3, 1) are predom- 
inant during the slow tail phase of the disturbance. It is that slow tail of the den- 
sity disturbance beginning about one hour after the onset of the storm v/hich 
contains most of the spectral energy and which is responsible for the observed 
world wide response of the thermospheric density during geomagnetic storms. 
Its dependence on storm time, latitude and longitude is discussed and compared 
with available satellite drag data. 
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THE RESPONSE OF THE THERMOSPHERIC DENSITY 
TO AURORAL HEATING DURING GEOMAGNETIC 

DISTURBANCES 


1. INTRODUCTION 

During geomagnetic disturbances the thermospheric density increases (Jacchia, 
1959), at least between 160 and 1000 km altitude (Jacchia et al. 1967). We shall 
refer to this phenomenon as "geomagnetic activity effect" of the thermospheric 
density. According to Jacchia and Slowey (1964a; 1964b), the intensity of the 
geomagnetic activity effect is proportional to the planetary index a^ for large 
disturbances and is proportional to the index for small disturbances. The 
effect occurs on a world wide basis. But the density disturbance appears to be 
systematically larger at high latitudes and during local midnight than at the 
equator during local noon (Roemer, 1970). The time lag between the peak of the 
geomagnetic disturbance and the corresponding density increase is about 5 to 7 
hours at low latitudes. That time lag is shorter by about one or two hours at 
high latitudes than at the equator (Jacchia et al, 1967). The density disturbance 
has an impulse form with a typical pulse width of one day. 

Thomas and Cning (1969) and Volland (1969) applying a simple one dimensional 
vertical model reproduced the height profile and the mean time lag of the density 
disturbance and concluded that the heat input responsible for the geomagnetic 
activity effect occurs predominantly within the lower thermospheric between 100 
and 200 km altitude and that the time delay is the natural response of the 
thermosphere to a pulse type disturbance of the heat input. 
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From the observations as well as from theory it is obivous that the most likely 
heating mechanisms of the geomagnetic activity effect are due to the penetration 
of high energetic charged particles into the aurora ovals thus ionizing the neutral 
gas within a small band near :i-05 geographic latitude, and energy transfer into 
heat of the neutral thermosphere via Joule heating (Cole, 1970) and recombination 
processes. These heating mechanisms are of course rather local phenomena. 
Therefore, the questions arise how the heat is transfered from the auroral regions 
to lower latitudes and why the response of the thermospheric density is worldwide 
in spite of the local heat input. We shall attempt to answer these questions in 
this paper with the help of a three dimensional thermospheric model, 

2. THE DISTRIBUTION OF THE HEAT INP\JT 

During periods of enhanced solar activity, fast ionized particles can penetrate from 
the magnetosphere into the aurora ovals , predominantly on the night time hemi- 
sphere (Akasofu, 1968). Tney react with the neutral gas mainly in the height range 
between 100 and 200 km and cause auroras and polar electro jets. We expect that 
the corresponding heating of the thermosphere occurs just in those ranges of height 
and geographic location during geomagnetic storms. 

Since the problem we are concerned here with is rather complicated we want to 
select a model which is as simple as possible. We therefore assume that the 
heating process is limited to a narrow band at ±65° geographic latitude on the 
night time hemisphere during equinox and we neglect the inclination of the 
geomagnetic dipole with respect to the geographic axis: 
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Here is A (in erg/cm ^ sec) a constant heat input at an altitude of ■' 100 
km within the strip between and h :,<(> latitude, 7 is the local time (zero at 
local midnight). For conveniance we assume very simple structures for the 
height distribution and the time variation of the disturbance: 
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where z is the height above ground and H is a constant scale height, and 

'oxp (-at) t ^ 0 

f (t) . for 

^0 t < 0, 

where t is the storm time and is a constant decay factor. 


(3) 


As it is well known (Chapman and Bartels, 1951) each regular function or a sphere 
can be represented by a sum of spherical surface harmonics. We develop the 
space distribution of Equ. (1) into a series of spherical harmonics: 
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is the local time which is related to a "zero" longitude, „ and is the midnight 
longitude during the onset of the geomagnetic storm, t “ 0, while \ is the usual 
geographic longitude with respect to the Greenwhich merdian. w is the angular 
frequency of the earth’s rotation. 0 is the total height and time integrated heat 
input ( in erg) - r,, f z ,,, where R is the earth’s radius and is the lower 
lx)undary of the heat input. 

Applying the well known properties of spherical functions (e. g. , Chapman and 
Bartels, 1951), we derive from Equs. (4) and (1) the coefficients 
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the total area where the heating occurs. 
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n, m 

0 , 0 

2, 0 

4, 0 

6, 0 

8, 0 

1, 1 

3, 1 

5, 1 

3, 3 

5, 3 

Qm 

1.00 

3.66 

2.22 

-2.65 

-6.91 

0.807 

3.58 

5.26 

-0.089 

-0.589 


Table 1: Spherical coefficient Qm of tlic heat input. For details see tejd;. 

Table I contains the numberical values of of the lower domains (n, m) for a 
co-latitude of = 25‘’ (9 q = 65°). We notice immediately that this series converges 
very slowly and that we need a large number of terms to approximate suffi- 
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ciently well our assumed heat input distribution of Equ. <1). Tnis number is 
determined by the broadness of the disturbed region and is of the order of n 
100 if the disturbed region has a thickness of about 500 km ( / 5”). We shall 

see in the next secticm, however, that the thermosphere "filters out" the energy 
modes with large wave domain numbers n such that only the energy modes with 
low wave domain numbers become significant generators of the geomagnetic 
activity effect on the thermospheric density. 

3. THE THERMOSPHERIC SYSTEM TRANSFER FUNCTION 

In this section we want to determine the response of the thermospheric density 
to a heat input of the general form of Equ. (1). In order to do such calculations, 
we remember from system theory that the output response of a system, f) (t), can 
be calculated from the product of the system transfer function, G { ,), with the 
input frequency spectrum Q (a): 

/-(t) ^ 6 ( 0 oi-t (1 V ( 5 ) 

*'•00 

where Q (w) is the Fourier transform of the input Q (t) which is in our special 
case the heating function of Equ. (4) (e. g., Stein and Jones, 1967). 

We can easily determine the frequency spectrum of the heat input of Equ. (4) 
which is in complex representation 
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It remains to determine the system transfer function G ( ) of the thermospheric 
density. In order to do this, we again Imve to restrict ourselves to a highly 
idealistic model. For conveniance, w'e shall assume that perturbation theory 
is valid. Under this condition, full wave theory at thermospheric heights 
has shown that harmonic density waves behave like quasi-evanescent waves 
at altitudes above about 200 km. Their vertical dependence therefore can be 
approximated by 


where H^, is a real scale height which is slightly greater than the scale 
height of the quiet thermosphere Ho (Volland, 1970). That means that these 
waves are attenuated very quickly outside of the location of thier generation 
and that one can neglect in a first approximation their propagation properties. 
Futhermore, we assume a time constant ion-neutral collision number 
of the order 


' 4 X 10 ’^ sec“i ^9) 

at 350 km altitude. This collision number is much greater than the Coriolis 
term 2 \l cos . Thus, we shall neglect the Coriolis force in the equations of 
horizontal momentum. 


We set 
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where p"*, , u™ and v™ are the coefficients in a series of spherical harmonics 

of pressure, density and horizontal winds of the wave field in the frequency 
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domain. . is the coefficient of molecular viscosity. In that approach the vertical 
gradients of the horizontal winds tend to become zero with increasing height 
which is reasonable. Furthermore, we assume isothermal conditions for 


the quiet thermosphere with the mean values of pressure and density 






( 11 ) 




and set 



Then, from the equations of horizontal momentum conservation it follows for 


each term (n, m) of the wave in the frequency domain at the angular frequency : 
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r = Rq + z (Rq = radius of the earth). 
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In nur approximation, the viscosity force behaves like a collision term. It is of 
the order 

^ Vi* *' 1 src”‘ 

at 350 km altitude and remains constant with height. In a more exact treatment, 
it turns out that . slightly increases with height and becomes the dominant 

V' I * 

friction term abfwe about 400 km (Geisler, 1967; Kohl and King, 1967). 


Using Equ. (13), we determine the horizontal divergence of the wind field as 
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Proceeding in a manner outlined in (Volland, 1970), we find from the equations 
of continuity and of the energy balance the following relationship between the 
frequency domains of heat input and density: 
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^0 3 x 10~3 

2 n (n 1) 


soc^/cm^] 


a Hj 'Hq 1.15 

'/ 1.5(ratio between the specific heats) 
c r 840 m/sec (velocity of sound) 

/< r. (coef f icients of the heat conduction). 


The numbers of R. and a have been adopted from full wave calculations done for 
the domain numbers (0, 0), (2, 0) and (1, 1) at a frequency of = for a 
thermospheric model of exospheric temperature of = 1000°K and assuming 
H = H 0 . They are valid within the height range between 300 and 400 km. These 
numbers will be kept constant, independent on n or u in our approach, though 
this is certainly a rough approximation. As we shall see in the following, 
the wave coefficients of the three above mentioned wave domain numbers play the 
predominant role for the geomagnetic activity effect. Therefore, we expect that 
our approximation gives rise to errors which are not too serious. 

The physical meaning of the different terms in l^qu. (15) is the following: 

Equ. (15) which is essentially the first law of thermodynamics contains two 
reversible terms which describe the adiabatic increase of internal energy of the 
gas (the term C) and the work done by the gas (the term L^), the last one mainly 
due to horizontal winds. The dissipative terms are due to vertical heat 
conduction (the term R.) and due to horizontal ion drag and viscosity (the term 
R^). Both terms behave like heat sinks which take energy out of the wave and 
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convert it into internal energy of the surrounding gas. It is apparent from Equ. 

(15) that at very low frequencies only the dissipative terms are effective. In the 
case of the zero component (n = 0) wave dissipation is solely due to vertical heat 
conduction. For the higher wave domain numbers (n > 0) horizontal ion drag and 
viscosity become increasingly important as dissipation mechanism because 
the horizontal scale length of the wave structui^e decreases with increasing n. 

The system transfer function {a ) of Equ. (15) which only depends on the wave 
domain number n has the form of the amplification factor of an R-amplifier with 
L-equalizer (the block diagram of which is shown in Fig. 1): 

(.) :: (16) 

^ u 

R 

A 

Here, u^ - QJJ R. is the grid voltage of a triode, R. is the internal resistance of 
the triode, u^ = is the anode voltage and R^, and C are resistor, inductor 
and capacitor of the anode circuit. The properties of such an amplifer are described 
in every textbook of system theory. 

In Fig. 2 we plotted the magnitude of versus frequency co with the wave domain 
number n as parameter. We notice that for small numbers n the transfer G 

n 

behaves like a low pass filter with bandwidth of 

^0 ~ 2.5 x 10"'* sec“^. (17) 

The anode resistor R and the inductor L decrease with n. Thus, the anode 
circuit becomes a short circuit with increasing n, and the amplification decrease 
like l/n^ for n > > 1. However, in the vicinity of 
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there occurs a resonance peak due to the over equalizing effect of the inductor 
L ^ . It gives rise to a maximum of of 
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That resonance effect can already be seen from a simple theory of plane internal 
gravity waves. In order to compare our result with plane wave theory, we notice 
that the zonal spherical functions PJ have n zeros between 0 n < 180°. Thus, 
their ’’horizontal wave length” is 


(n > > 1) 

X n pj 


and their normalized ’’horizontal wave number” for the resonance frequencies 
Q is 
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The range 1 < S < 2 is just the range where the attenuation of plane inter. .al 
gravity waves has it minimum (e. g., Volland, 1969b). However, its deepness 
is much smaller at altitudes above about 300 km than Fig. 2 suggests and in fact 
nearly disappears at 400 km altitude. Therefore, we must conclude that our 
approach, though being qualitatively consistent with the exact theory, overesti- 
mates the amplification effect of the thermospheric circuit for great wave domain 
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numbers n. The reason for that discrepancy lies in our assumption of constant 
circuit elements R. and n (n + 1) , on the one hand. One the other hand, 

viscosity and heat conduction due to horizontal gradients of winds and temperature 
which are neglected in our approach lead to an additional decrease of for great 
domain numbers n. Our whole concept of course breaks down for accoustic wave 
( - > 10" 2 sec"^ ). 


The resonance frequencies lie well beyond the collision frequency i . 
Therefore, the error due to our approach involves predominantly the high fre- 
quency band of the disturbance. As we shall see in the next section, the 
decay factor a in Equ. (3) related to an average geomagnetic disturbance is 
of the order 

a ,Q < < .0^, (n > 3), 

and the bandwidth of the spectral function in Equ. (8) is therefore 

- Q. 

Thus, the low frequency tail of the disturbances which contains most of the spec- 
tral energy is not affected very much by the high frequency range . 

Because of these reasons, is appropriate to make a further simplification by 
taking the inductor = 0. That gives the simplified system transfer function 
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Now, ( *;) resembles the reziprocal impedance of a simple RC-circuit. 

The magnitudes of are plotted as dashed lines in Fig. 2 and indicate that only 
the high frequency band of G^ for great n is changed by that simplification. 

4. THE DENSITY RESPONSE AS FUNCTION OF TIME 

In this section we want to determine the density disturbance in the time domain 
caused by the heat input of Equ. (4) and by the system transfer function G('-) given by 
Equ. (21). As already mentioned, by using G (a) instead of G (o ) of Equ. (15) 
we suppress the high frequency response of the density disturbance and consider 
only the slow tail of the disturbance. The transformation from the frequency 
domain into the time domain using Equ. (6) leads to a density distrubance of 


with 


A 


a Q 

4 rr rj H 


g ('z') 



/C Pn > 


(22) 


Qm 

An (A, t) = — H_ {e““* cos (tn r - a™) 
A™C 


e * cos (tn (k - k^) - a™)} 


and 


A™ = - a) 2 + m2 n2 



Since A” increases with n and m , the coefficients with large domain numbers 
(n, m) are filtered out rather rapidly within the thermospheric "circuit system". 
The density distribution of Equ. (22) not only depends on the storm time t but also 
on local time ^ . We shall discuss first that part of the density disturbance which 
only depends on storm time (m ~ 0). 


4,1 The Zonal Components 

For the zonal components (m = 0) it is 
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That, time of the maximum shifts toward lower values as x^, increases. From the 
numerical values in Equ. (15) it follows 


= 1/(R. C) ~ 2.5 X 10“5 sec"i. 
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From the observed averaged time lag between geomagnetic disturbances and 
the thermospheric geomagnetic activity effect of hours we determine 


and from Eq. (24) 


^0 ^^mnx^O 
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The time of maximum of the/. -coefficient then becomes 


^^max^2 ^300 sec"* ' l.l hours, 


and it is 


('t < <'!• 

' mnx'n ' mox'n“l* 

Thus, the wave components with wave numbers n > 2 are significant only within 
the first hour after the commencement of the disturbance and then decay rapidly. 


In Fig. 3a we plotted the coefficients (t) of the density versus storm time t 
for the wave domain numbers 0, 2, 4, 6, and 8 using the values of Q° from Table 1. 
and the number, a, from Equ. (25)*. We note that in spite of the larger coefficients 
Q® as compared with QO = 1, the zero coefficient po (t) dominates about 3 hours 
after the beginning of the disturbance. 

The two first coefficients p® and p® already dominate 30 minutes after the beginning 

V A 

of the disturbance. However, the higher harmonics with n > 2 may produce signif- 
icant amplitudes within the first 30 minutes. Therefore, we expect an oscillating 
behavior at the beginning of the disturbance which can not be produced in our 
simlified theory. 

* Positive values are indicated by solid lines, negative values are indicated by dashed lines. 
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The sum of the zonal components is plotted versus storm time t in Fig. 3b for 
four different latitudes % =0°, 30", GO" and 90' . Here we note a shift of the 
maximum disturbance in the density from 6 h at the equator (cji ® 0°) to 2,5 h at 
the poles (9 = ± 90°) and an increase of the maximum density amplitude by a 
factor of about two between equator and poles. Fig, 4 presents that latitudinal 
dependence of the zonal maximum amplitudes and the zonal maximum times versus 
latitude. That latitudinal behavior is essentially due to the influence of the zonal 
harmonic which is positive at latitudes | (p| > 35° and negative at lower latitudes. 

The total duration of the impulse is of the order of two days. At the equator the 
initial phase is negative during the first hour. Ten hours after the start of the 
disturbance the amplitudes of the zonal components with n > 0 are decayed 
and there remains only the zero components / q Pg . 

In order to study the influence of the shape of the heat input on the density 
disturbance we varied the decay factor a of the heat input and calculated the 
corresponding times of maximum and the maximum amplitudes 

of the zero component. We plotted a and (Pg)^a^ versus (t^^^ )o in Fig. 5a. 
Moreover, we determined in Fig. 5b the total duration AT of the impulse versus 
maximum time. That total duration is not uniquely defined because of the 
exponential decay of the density disturbance. Therefore, we define that the end 
of the disturbance is reached where the amplitude has dropped to 

pg (t :: A T) ■- 3 X 10"® secVcm^ 

In our example in Fig. 3b that condition is equivalent to a ratio 
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That ratio becomes greater with decreasing a-fector and vice versa (see Fig. 

5a) and shall take account of the fact that in practical data processing the 
evaluation of a disturbance is limited by an absolute noise level rather than 
by a relative amplitude ratio. Fig. 5s shows that the total duration ^ T 
increases with increasing time maximum (t^„^)5 . If we use any other defi- 
nition for the duration of the disturbance we change the slop of the curve in Fig. 

5b without changing the trend. 

The results of Fig. 4 and 5 should be compared with available satellite drag data. 

The decrease of the time of maximum with increasing latitude is in agreement 

with the analysis of Jacchia et al., (1967) and Hoemer (1967). The increase of 

the maximum amplitude by a factor of about two between equator and poles agrees 

with Roemers (1970) analysis. However, in Roemers (1970) paper, he did not 

confirm the latitudinal dependence of the mean values of the maximum time but 

claims that there is no significant dependence* . We shall return to that discrep- 
# 

ancy between observation and theory in section 4.2. Roemer (1970) also deter- 
mined a relationship between duration of the impulse and the time of maximum 
and found a linear dependence (the dash-dotted line in Fig. 5b). That observation 
is sufficiently well reproduced by our theroy. 

Fig. 5a also shows a relation between the decay factor p and the maximum 
amplitude maximum decreases rapidly with increasing, a, indicating 

that isolated short periodic polar substorms lasting typically one hour (a ~ 10~^ 
sec"^) do not generate a significant geomagnetic activity effect. 

* Though his individual data scatter widely from 0 to 15 hours within the whole range of latitude. 
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We can estimate the amount of heat input necessary to cause the observed 
density disturbance during magnetic storms. From Hoemers (1970) analysis, 
we determine an Increase of the exospJieric temperature of ‘ T, - 50 K at a 
latitude of i -- 0 for a geomagnetic index of * 3. That increase of the 
exospheric termperature in turn is related to u maximum relative density 
amplitude of 


or 


■ *■ trax^' 0 0,17 


max IS ' 


(27) 


at moderate solar activity and at 3f-0 Ian altitude. The total heat input integrated 
over time and space is then according to Equ. (22) 
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4 " r2 H 
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K (z) n 



max 


where r^ = 0500 Ian, z = 350 km and H = 50. From the solid curve with 
parameter = 0^' in Fig. 4, a maximum value of 


1.4 X lO"*^ sor3/em2 

n n I 

J max 

was used. 

The total heat input of Equ. (28) s-hould be compared with an estimate made by 
Akasofu (1908) for the energy content carried by the auroral electrons and the 
ring current protons during a polar substorm. His estimate is 2 x 10^^ erg. 
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Comparing this number with our value in Equ. (28) we have to conclude that a 
moderate geomagnetic storm should contain at least three times the energy of 
a substorm if we assume an efficiency factor for the heat transfer of • < 0.3. 

Within the auroral belts the maximum height integrated heat input is according 
to Eq, (5a) 


H ti Oq Q p ~ 20 erg ^(ctn^ sec) 


where the numbers F “ 1 x 10 cm*; « 25® and A * 5® have been adopted. 
That value seems to be not unrealistically large. The maximum heat input at 
t * 0 within a vertical column averaged over the whole sphere is 


a 9 / (4 77 f2) '' 0.4 c rg/(cm* sec)« 

That value is only a fraction of the XUV-input into the thermosphere above 100 
km (Hinteregger et al. , 1965). 


4.2 The Local Time Dependence of the Density Disturbance 

We now turn to the discussion of that part of the density disturbance in Equ. (22) 
which depends on local time r . This part is connected with the associated 
spherical functions P™ (m > 0). In Fig. 6 we plotted the coefficients (m > 0) 
versus storm time from Equ. (22) adopting the number a from Equ. (25) and the 
coefficients Qjj from Tab. I. Fig. 6a is calculated for the longitude ofX-\Q=0°, 
and Fig 6a is valid for a longitude \ - 90°. Tne curves for the longitude of 

180° and of 270° differ from those of Fig. 6a and 6b, respectively, by only a 
change in sign. Positive values in Fig. 6 are indicated by solid lines, negative 
values by dashed lines. From Fig. 6 we note that the coefficients with larger 
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wave numbers n peak earlier than those with lower wave numbers, a fact wliich 
we observed have from the belmvior of zonal coefficients in Fig. 3, Moreover, 
the coefficients with m = 3 are of no signifiance throughout the temporal range 
considered. This again is due to the small excitation factors of waves with 
large wave domain numbers (n, m). Two hours after the onset of the distrubance 
the predominant coefficient remains / |. During its maximum phase that coeffi- 
cient is positive at local night and negative at local noon. Thup it contributes to 
the total density disturbance in such a manner that it increases the effect during 
the night and decreases it during the day. This behavior becomes more clear 
in Fig. 7 where we plotted the total density disturbance X versus storm 

time for the three latitudes cp = 0^, 30^’ and 60"^ and for four different longitudes 
,V - \ " Of, 90°, 180° and 270°. In Fig. 7 we omitted the curves valid at the 
poles (/' = 00°) because they are identical with the corresponding curve in Fig. 
3b and show no dependence on local time. 

As is already obvious from Eq. (22), the shape of the disturbance, its maximum 
amplitude and maximum time depend on storm time, local time and individual 
longitude of each event. 

The maximum amplitudes and the times of maximum for the total density dis- 
turbance are plotted versus longitude in Fig. 8a and 8b for the four latitudes 
considered. The maximum amplitude shows a quasi-harmonic variation v/ith 
longitude with maximum values about \ - a.^ = 330°. The amplitudes of the 
curves in Fig. 8a are largest at 60° latitude, the amplitudes of the maximum 
times in Fig. 8b are largest at the equator. If we had plotted these curves versus 
local time we had to shift the abszissa in Fig. 8 by H t^^^^ . Then, the maximum 
of the curves in Fig. 8a occurs at about 
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•w f 330'’ - O ' to 90“ (Oh to 6h Lt) 

that is during the early morning hours, if we consider the range of the times of 
maximum from Fig. 8b to be 2h to 8h Lt, 

We compare the result of Fig. 8 with observations of Roemer (1970). Roemer 
found a dependence of the maximum amplitudes on local time at low latitudes 
which is consistent in magnitude and phase with our result of Fig. 8a. However, 
he did not observe a signicant longitudinal dependence of the maximum times as 
we predict in Fig. 8b. A possible explanation for that discrepancy as well as for 
the similar discrepancy concerning the latitudinal dependence of the maximum 
times (see section 4.1) may be the fact that each individual geomagnetic activity 
effect depends on its own midnight longitude Observations at longitudes 
arbitrarily related to that a ^ lead to a wide scattering of data for the maximum 
times between 1.4 and 8 hours in our special case of Fig. 8b. Moreover, these 
data depend on the decay factor, a, of every individual storm. Finally, the prob- 
lem becomes even more complicated by the facts that (a), the temporal behavior cf 
any individual storm may be far from similar to our assumed simple exponential 
model of Eq. (3) and (b), that with increasing intensity of the storm the auroral 
ovals shift to lower latitudes , thus influencing the ratios between the spherical 
coefficients of the heating function in Eq. (4). At present, therefore, these dis- 
crepancies remain open questions. We would suggest that geomagnetic activity 
should be studied in terms of the midnight longitude of each individual storm. 

5. DISCUSSION AND CONCLUSION 

In the foregoing sections we considered on impulse type heat input within a small 
band of latitude near the auroral ovals during local night. We determined the 
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"system transfer function" of the thermospheric density and calculated from that 
function and from the Fourier spectrum of the heat input the temporal response 
of the thermospheric density to the heat input. It appears tliat the components 
with large wave domain numbers (n, m) have excitation factors which decrease 
rapidly with increasing numbers (n, m). Thus, after the onset of the 
geomagtietic storm, the density distribution of the geomagnetic activity 
effect can be described by the four spherical harmonics with the lowest wave 
domain numbers. Three hours after the onset of the disturbance the zero compo- 
nent (0, 0) dominates and leads to a slow tail of the density amplitude and thus to 
the observed global response of the thermospheric density. We examined its 
dependence of this effect on storm time, latitude and longitude and compared the 
theory with available satellite drag data. While the dependence of the maximum 
amplitude on latitude and longitude could be well reproduced in this theory, a 
discrepancy remains between observations by Roemer (1970) and our calculations. 
This concerns the dependence on latitude and longitude of the times of maximum 
for geomagnetic activity effect which is not significant according to Roemer. We 
suggest that this discrepancy may be due to the dependence of the geomagnetic 
activity effect on storm time which has not been taken into account in an appro- 
priate manner in previous treatments. 

Our theory does r / appropriately describe the high frequency range of the 
geomagnetic acti\'"y effect which might be significant during the first hour of 
the disturbance, Thus we have excluded short periodic gravity waves and 
accoustic waves which have in fact been observed in the vicinty of the auroral 
zones just after the beginning of the geomagnetic disturbance (Taeusch and Carignan, 
1970; Tastud, 1970). Tnese disturbances can not be detected from satellite drag 
data due to the poor time resolution of this method which is not better than 
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three hours. Also, the short periodic geomagnetic bay disturbances or polar 
substorms do not supply enough spectral energy into the slow tail to become 
detectable by the satellite drag method. Those short period disturbances may 
be described more appropriately in terms of plane gravity wave propagation 
(e. g. Chimonas and Hines, 1970). 

Our theory is appropriate essentially for the explanation of the satellite drag 
data which observe the slow tail of the geomagnetic activity effect. The slow 
tail however contains most of the spectral energy of an event with a typical 
duration of one day. 

The local time dependent component (1, 1) of that slow tail gives rise to horizontal 
winds which are similar in structure to the winds of the tidal (1, -1) mode at 
lower altitudes which are however shifted in phase by about 120“ with respect 
to the tidal wind. In the light of that result one should reconsider the orgin of the 
geomagnetic SD-current in terms of djmamo theory. As Fulmshima and Oguti (1953) 
showed, dynamo theory can explain the SD-current if the wind is shifted in phase 
by 100“ to 150° with respect to the wind responsible for the Sq-current. Since 
the wind of the Sq-current is the wind of the tidal (1, -1) mode (Stening, 1969) 
and since the phase of the wind of the geomagnetic activity effect is shifted in the 
right manner, we suggest that at least part of the geomagnetic SD-current may be 
generated within the thermosphere by the dynamo action of winds related to 
the geomagnetic disturbance. 
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FIGURE CAPTION 


Fig. 1 Block diagram of an R-amplifier with L-equalizer 

Fig. 2 Magnitude | G„ ! of the "signal transfer function" of the thermospheric 
density versus angular frequency . , The wave domain number n is 
shovTi for each curve. The dashed curves give the approximate function 
I G|J . For details see text. 

Fig. 3a Zonal eoefficients / ® versus storm time t. Full linos indicate positive 
values, dashed lines indicated negative values. 

Fig. 3b Sum of zonal components / [J PJ] versus storm time t for four 
different latitudes <p . 

Fig. 4 Maximum amplitudes (solid lines) and times of maximum (dashed 

lines) of the sum of the zonal components (Fig 3b) versus latitudes cj.-. 

Fig. 5a Maximum amplitude (/ and decay factor a of the heat input versus 
maximum time (t^gx)o zero component. 

Fig. 5b Total duration A T of the density disturbance versus time of maximum 
(tmax)o zero component. 

Fig. 6 Spherical components /Ojj P™ (m > 0) versus storm time t calculated 
at the equator (9= 0°). Solid lines indicate positive values, 
dashed lines indicate negative values. 

a) Longitude \ - \ = 0*^ 

b) Longitude A. - = 90° 
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Fig, 7 Sum of all components ^ P” versus storm time t for three 
different latitudes 9 and four different longitudes q . 

Fig. 8 Maximum amplitudes (Fig. 8a) and maximum times (Fig 8b) taken 
from Fig, 7 versus longitude \ for four different latitudes 9. 
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TRANSFER FUNCTION lOn (sec3/cm2) 


ANGULAR FREQUENCY w (sec-‘) 

Figure 2. Magnitude iGnjof the "signal transfer function" of the thermospheric density versus 
angular frequency u-. The wave domain 'number n is shown for each curve. The dashed curves 
give the approximate functions iG„l. - For details see text. 
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Figure 3b. Sum of zonal components - . ^ P ^ versus stom time t for four different latitudes 
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